FERMIL AB-PUB-07- 105- A 



Impact of astrophysical processes on the gamma-ray background from dark matter 

annihilations 
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We study the impact of astrophysical processes on the gamma-ray background produced by the 
annihilation of dark matter particles in cosmological halos, with particular attention to the con- 
sequences of the formation of supermassive black holes. In scenarios where these objects form 
adiabatically from the accretion of matter on small seeds, dark matter is first compressed into very 
dense "spikes", then its density progressively decreases due to annihilations and scattering off of 
stellar cusps. With respect to previous analyses, based on non-evolving halos, the predicted anni- 
hilation signal is higher and significantly distorted at low energies, reflecting the large contribution 
to the total flux from unevolved spikes at high redshifts. The peculiar spectral feature arising from 
the specific redshift distribution of the signal, would discriminate the proposed scenario from more 
conventional astrophysical explanations. We discuss how this affects the prospects for detection and 
demonstrate that the gamma-ray background from DM annihilations might be detectable even in 
absence of a signal from the Galactic center. 

PACS numbers: 95.35.+d, 97.60.Lf, 98.62. Gq 



I. INTRODUCTION 

Indirect dark matter (DM) searches are based on the 
detection of secondary particles such as gamma-rays, 
neutrinos and anti-matter, produced by the annihilation 
of DM particles either directly, or through the fragmen- 
tation and/or decay of intermediate particles (for recent 
reviews see Refs. [J HIl). 

Among the proposed strategies of indirect detection, 
searching for a diffuse gamma-ray background produced 
by the annihilation of DM in all halos at all redshifts 
appears particularly interesting, because of the useful in- 
formation that such a signal would provide on the dis- 
tribution and evolution of dark matter halos [J, H, @- 
Previous calculations have been performed under the hy- 
pothesis that the shape of DM profiles doesn't change 
with time, a circumstance that led to the conclusion that 
the prospects of detecting gamma-rays from the Galac- 
tic center (GC) are more promising than the gamma-ray 
background [7j. However, the annihilation signal mainly 
comes from the innermost regions of the DM halos, i.e. 
regions where the gravitational potential is dominated by 
baryons, and where the extrapolation of numerical sim- 
ulations is most uncertain. 

In particular, the strong evidence for supermassive 
black holes (SMBHs) at the centers of galaxies suggests 
that the DM profile is inevitably affected by astrophysi- 
cal processes on scales that cannot be resolved by numer- 
ical simulations [8j . The formation of massive black holes 
(BHs) at the centers of DM halos can significantly modify 
the DM profile, especially if the process of BH formation 



happens "adiabatically", i.e. the formation timescale is 
much longer than the dynamical timescale of DM parti- 
cles around it d, El El, Ell ■ These so-called spikes 
of DM inevitably interact with stars and other structures 
in the Universe (e.g. binary black holes), a circumstance 
that typically leads to a decrease of the DM density, and 
thus of the annihilation signal [H, El, El, E3 • 

In order to detect the enhancement of annihilation ra- 
diation from these dense structures, one thus has to look 
either for spikes where astrophysical processes are less ef- 
fective, that evolve in regions with low baryonic densities, 
as in the case of intermediate-mass black holes [18j |. or 
for the contribution to the gamma-ray background from 
spikes at high redshift, when the DM enhancements had 
not yet been depleted by astrophysical processes. 

It is therefore important to re-analyze the prospects for 
detecting the gamma-ray background produced by cos- 
mological DM annihilations, in a self-consistent scenario 
that takes into account the time-dependent effect of as- 
trophysical processes on the distribution of DM. Here, we 
first provide a prescription to assign BH masses and stel- 
lar cusps to generic halos of any mass and at any redshift. 
We then follow the formation of spikes around SMBHs 
at high redshift, and their subsequent disruption due to 
the interaction with the stellar cusp and to DM angula- 
tions. Finally, we integrate the annihilation signal over 
all redshifts and all structures and discuss the prospects 
for detecting the induced gamma-ray background. 

The paper is organized as follows: In Sec.|TT]we specify 
how to assign spikes to cosmological halos of given mass 
and a given redshift, and how spikes evolve. In Sec. lIIII we 
calculate the gamma-ray background produced by DM 
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annihilations in halos of all masses and at all redshifts. 
Finally in Sec. lIVI we present our conclusions. We include 
the description of the halo density profile, its mass dis- 
tribution, and evolution for the sake of completeness and 
to allow comparison with existing literature in Sec. 
Throughout this paper, we assume a flat ACDM cosmol- 
ogy with Q m = 0.3, h — 0.65, spectral index n = 1, and 
er 8 = 0.9. 



II. ASSIGNING SPIKES TO HALOS 

To estimate the effect of BHs on the gamma-ray back- 
ground produced by DM annihilations, we first need to 
model the formation and evolution of BHs in halos of 
given mass and at a given redshift, and to follow the 
formation of DM spikes, and their subsequent destruc- 
tion due to scattering off stars and to DM annihilations. 
Strong constraints on the BH population at all redshifts 
come from the relationships between DM halo proper- 
ties and BH masses observed in the local universe, and 
from the quasar luminosity function. In this section we 
devise a strategy to assign BH masses to host halos at 
any redshift, and to calculate the DM distribution in the 
resulting spikes. Further details on the normalization of 
DM halos, and on the their cosmological evolution, can 
be found in the Appendix. 



A. SMBH formation 

In ACDM cosmologies, DM halos (~ 10 8 M ) begin to 
form at large redshifts (z ~ 20) and subsequently grow 
through mergers, while stars form from gas that falls into 
the halo potential wells. At some point, SMBHs form 
from the stars and gas at the centers of the halos. Exactly 
how this occurs is not clear. However, the luminosity 
function of quasars as a function of redshift traces the 
accretion history of these BHs [l9| , suggesting that BHs 
grew significantly, by accretion, from their initial seeds, 
with large mass-to-energy conversion efficiency [2(i l2~ll . 
I22I ]. An estimate of the average growth history of BHs 
presented in Ref. [12] , suggests that the redshift by which 
BHs have reached 50% of their current mass, varies with 
the BH mass, ranging from z > 2, for BHs more massive 
than 10 10 Af o , to z < 1 for BH masses below 1O 6 M . 
We adopt here a simplified approach, where all BHs were 
already in place at a characteristic redshift of formation 
z = zbh, and will discuss the dependence of our results 
on z B H- 

In the local universe, tight empirical relations are ob- 
served between SMBH mass and the mass of the DM 
halo [23| and the luminosity [24[ and velocity dispersion 
[IH of the stellar component. Based on these results, we 
adopted the following prescription for assigning SMBHs 
to halos: 

1. The local correlation between SMBH and halo mass 
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FIG. 1: Msmbh as a function of the host halo mass Matz = 
6 (dashed lines) and z — (solid lines). The three relations 
in Eq. (fTJ) (a)-(c) are shown from bottom to top for both 
redshifts. 



[231 ] is used to calculate the mass of the SMBH 
(Msmbh) lying in a halo of mass M at z = 0. 

2. A SMBH with this mass is placed in the progenitor 
of this halo at z = zbh ■ 

3. The halo is evolved from redshift zbh to 26], 
while leaving Msmbh fixed. 

Based on Ref. (23|, we considered the following rela- 
tions between Msmbh and M at z = 0: 



Msmbh 
10 8 M n 



0.027 (M^o) 1 - 82 
0.10 (M^o) 1 ' 65 
0.67 (Mho) 1 '" 



(a) 
(b) 
(c) 



(1) 



where M 12 ,o = M(z = 0)/10 12 M Q . The differences 
reflect different assumptions between the virial radius 
r v and the circular velocity. Figure [T] shows the above 
three relations between Msmbh and M. The Msmbh 
obtained at z — is subsequently placed in halos at 
z = zbh- For a fixed halo mass, SMBH at z = is 
less massive than SMBH at z = 6. This is because the 
halos at z — 6 would have evolved to a more massive halo 
by z = 0, where the Msmbh is determined. 



B. Formation and evolution of DM Spikes 

The growth of SMBHs inevitably affects the surround- 
ing distribution of DM. In fact, it can be shown that the 
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adiabatic growth of a massive object at the center of a 
power-law distribution of matter with index 7c , induces 
a redistribution of matter into a n ew, steeper, power-law 
with index j sp = 2 + 1/(4 - 7c ) @, Hfl, [U Gj, . Such a 
DM enhancement is usually referred to as a "spike" . For 
the widely adopted Navarro, Frenk and White (NFW) 
profile (see Appendix for further comments and refer- 
ences), 7 C = 1, and the spike profile, immediately after 
its formation (i.e. at t = tf, the time when the spike is 
formed), can be expressed as 



Ps P (r,0) = p{r b ,o) 

Jb,0 



-7/3 



(2) 



inside a region of size r^o w 0.2 [271, where is the 
radius of gravitational influence of the SMBH that is de- 
fined as 



GMsMBH 



(3) 



where G is Newton's constant and a the one-dimensional 
velocity dispersion. Msm bh can be related to a through 
the empirical relation [25| 



MsMBH 

10 8 M m 



(1.66 ±0.24) 



200 km s" 



4.86±0.43 



(4) 



Eq. ([4J is known to be valid for SMBHs in the mass 
range 1O 65 M < M S mbh < 10 9 5 M o and may extend 
to higher and lower masses |25[ . 

Once the spike is formed, several particle physics and 
astrophysical effects tend to destroy it (e.g. 0, [IH, [H, 
\v\). Here we focus on the gravitational interaction be- 
tween DM and stars near the SMBH, which causes a 
damping of the spike, and on self-annihilation of DM near 
the SMBH, which decreases the maximum density of the 
spike. 

The DM and baryons gravitationally interact with each 
other. Stars in galactic nuclei have much larger kinetic 
energies than DM particles, and gravitational encounters 
between the two populations tend to drive them toward 
mutual equipartition. DM is thus heated up, dampening 
the spike while maintaing roughly the same shape of the 
density profile. Based on the results in Refs. [l|| Il6j |. 
we adopted the following approximate expression for the 
decay of the spike intensity with time: 



p(r, t) w p(r, 0) k ; k 



(5) 



where r is the time since spike formation in units of the 
heating time T heat [l]| 



T heat = 1.25Gyrx 

MsMBH 

3 x 10 6 M Q 



(UlY 

V 2pc / I ™* 



15 
hiA 



(6) 



in* is the effective stellar mass, and is equal to ~ 1.8Af 
assuming a Salpeter mass function and O.O8Af < < 




10-' 



FIG. 2: Redshift evolution of the spike parameters 
(solid lines) and X), t o (dotted lines), for various halo masses: 
Mh(zbh) = 10 13 , l0 12 , 10 11 M Q from top to bottom (red, 
blue, green), and assuming zbh = 6. 



2OM . In A = ln(0.4A), with N « 6 x 10 6 the number 
of stars within for our Milky Way Galaxy. Although 
r/i is a function of MsmbHi we approximate In A to be 
constant as Theat dependence on N is logarithmic. 

The size of the spike decreases, with respect to the 
initial value r b o > with time 



r b (t) = K S r bfil 5 = (7 sp -7c) 1 
The spike density profile is thus given as 



(7) 



>(r,t) 



.,0 K 



-lsp 



Ic/ilc-lsp) ■ (8) 



Figure [2] shows the spike evolution under the effect of 
DM interaction with baryons, assuming zbh = 6 for dif- 
ferent halo masses. The spike parameters are shown in 
units of the halo reference radius r$ (defined in Sec. [A"]) , 
so that x b = r b /ro and x bi o = r b ,o/rQ. It can be seen 
that spikes formed in very massive halos are not affected 
by heating from the baryons, whereas those in less mas- 
sive halos quickly dissipate. Hence low mass halos give 
negligible contribution to the gamma-ray signal. 

A robust lower limit on the size of the spike is provided 
by the last stable orbit (ri so ) of a test particle around 
the SMBH. However, annihilation itself sets an upper 
limit on the DM density. The evolution equation of DM 
particles at radius r and time t is 



(av) n sp (r, tf 



(9) 
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where the dot denotes a time derivative. Although this 
expression is correct for circular orbits, a more sophisti- 
cated approach would take into account the eccentricities 
of orbits, and would start from the single-particle dis- 
tribution function f(E,L) describing the DM particles, 
where E and L are the energy and angular momentum 
per unit mass respectively, and compute orbit-averaged 
annihilation rates. Such a calculation has apparently 
never been carried out and is beyond the scope of this 
paper. Under the assumption of circular orbits, one finds 
that the maximum number density at a given time t can 
be expressed as 



,{r, tj 



(10) 



1 + n sp (r,t f ) (t - t f ) (av) 
This is usually simplified to obtain a maximum density 

m v 1 



P P l(t) 



(av) (t-t f ) 



(11) 



The radius where p sp reaches this value, denoted as r p , 
can be calculated by inserting Eq. ([8]) into the above 
equation. The maximum allowed density decreases with 
time due to self-annihilation and a plateau of constant 
density forms from r p down to ri so . As r m in is larger 
than riso, except for very massive halos, the fully evolving 
spike density profile is given as 



Ps P (r, t) 




(riso < r < r p ) 
e(r\-^» , „ x -( 12 ) 

Figure [3] shows the density profile of an evolving spike 
which has formed at zbh — 6. Note that the evolution 
of the halo itself has not been taken into account and the 
halo and SMBH mass are fixed to 10 12 M Q and 10 7 M , 
respectively, at all redshifts in order to show only the 
changes due to the evolving spike. The halo profile is 
plotted in dashed line and the spike profiles at various rs 
are plotted in solid lines. The DM profile is divided into 
three regions; a plateau with magnitude p p i from r; s0 to 
r p , the prominent spike that scales as r ~ r ' sp from r p 
to rb, and the prominent halo with r ~ r _7c from to 
r v . Numerical computations, such as Ref. [l6| which has 
been calculated for our Galaxy, show the same features 
but with a smoother transition at rj,. 

It is convenient to express p sp ,o m terms of ro and 
Po, thus of halo mass. Given the halo mass M, one can 
obtain r v from Eq. ([AT]) and r fr om the definition of c 
of Eq. (|A6[) . It should be noted that ro is not a constant 
but varies with z and M. Similarly, po, which is also 
dependent on z and M , is obtained by solving 



M = 4-7T 



drr 2 p h (r) . 



(13) 



The reference spike density p sp ,o is normalized by the 
halo density at r = r&.o and r = 0. For a NFW halo this 
gives 



Po 



x hfi (i + x bt0 y 



(14) 




r (Mpc) 

FIG. 3: Density profile of an evolving spike. Halo and SMBH 
mass are fixed at 10 12 Mq and 10 7 Mq , respectively, at all 
redshifts, and zbh = 6. The dashed line is the halo profile, 
and solid lines are the spike profile in various stages of evo- 
lution. Spikes are plotted at r = 1, 2, 3, and ~ 14 (z = 0), 
from right to left (in order of decreasing size). The spike 
becomes prominent at rt(t), which decreases with time due 
to DM-baryon scattering. Self-annihilation of DM causes the 
maximum spike density value of p v i to decrease with time. A 
constant density of DM with value p p i is maintained from r p 
to the last stable orbit (~ 10" 12 Mpc). 



where Xb,o = r b.a/ r o, thus 



}(x) = Po 



Tap — 1 

C h,0 



(1 + Xbfi)'- 



KX 



(15) 



The expression of the total density profile depends on 
the redshift and radius. The total density profile is given 
as follows: 



Ptot 



Ph(r) (z > z B h) 

Ph(r) (z < zbh, r > r h ) 

Ph(r) + p sp (r,t) 

w p sp (r,t) (z < z B H, r < r b ) 



(16) 
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III. GAMMA-RAY BACKGROUND FROM DM 
ANNIHILATIONS 

The contribution of DM annihilations to the gamma- 
ray background flux <E> can be expressed as 0] 



HE) 



1 



4tt H, 



o 



a. 



where c here is the speed of light, Ho is the present value 
of the Hubble parameter, E s — E(l + z) is the energy 
emitted at the source, and h(z) = \/ (1 + z) 3 £l m + f^A- 
To allow an easy comparison with existing literature, we 
adopt a simple analytic ht to the continuum gamma-ray 
flux emitted per annihilation coming from hadronization 
and 7T° decay [I[ 



dN(E) 
~d~E 



0.42 



-8E/m, 



(E/m, 



11.5 



0.00014 ' 



(18) 



which is valid for E < m x . The exponential in the in- 
tegrand takes into account the effect of gamma-ray ab- 
sorption due to pair production on background photons. 
Following we write it as 



e~ T ^ « exp 



3.3(£/10GeV)- 



(19) 



The dimensionless flux multiplier £(z) can be written 
as the integral over all masses of an auxiliary function 
g(M,z), weighted by the halo mass function dn/dM, 
which is typically calculated in the framework of the 
Press-Schechter or Sheth-Tormen formalisms described 
in Sec. El 



CITI 

dM g(M, z) . 



(20) 



The auxiliary function g(M,z) is simply the flux multi- 
plier relative to a halo of mass M at redshift z, 



g(M,z) 



1 



dVp 2 {r), 



(21) 



normalized to the comoving background density squared. 
V is the halo virial volume, which is a function of redshift, 
and of the halo mass and concentration (see Eq. |A1[) ). 

The integration over DM spikes requires particular 
care. Since we are assuming that SMBHs do not evolve 
after their formation redshift zbh, the halo parameters 
in the £(z) calculation have to be evaluated at zbh, 
while the spike evolves with redshift as discussed above. 
Furthcmore, the M — Msmbh relationship must evi- 
dently break down at small masses. Here we have re- 
stricted the anlysis to spikes produced by BHs with mass 
M$mbh > 100 M Q , and have verified that the result is 
insensitive to this lower mass cutoff. 




FIG. 4: Gamma-ray background produced by DM annihi- 
lations in DM halos with spikes (solid lines), compared to 
the halo-only contribution (dotted). The EGRET diffuse flux 
limits [28| are shown for comparison. Halos have mass be- 
tween 10 s - 10 14 M e and are distributed according to the 
Press-Schechter formalism. The three expressions in Eq. (fl]) 
are used for the Msmbh-M relation: (a) (red) (b) (blue) (c) 
(green) from the bottom to top. The DM parameters adopted 
are m x = 100 GeV, (av) = 10 _26 cm 3 s _1 . 



Figure [4] shows the enhancement of the gamma-ray 
background due to the presence of spikes, compared with 
the standard calculation (halo only). For the figure, we 
have focused on the Press-Schechter formalism, but we 
find similar results for the case of ellipsoidal collapse a la 
Sheth & Tormen. We given an upper limit of z — 18 to 
Eq. (jTTJ) and assume that spikes form at Zbh = 2. DM 
parameters m x — 100 GeV, (av) = 10~ 26 cm 3 s _1 are 
used in the calculation. All three mass relations between 
the SMBH and halo of Eq. (fTJ) are shown in the figure, 
(a) to (c) from bottom to top in solid lines. The diffuse 
EGRET flux [28[ is plotted as a comparison. Enhance- 
ment due to the presence of evolving spikes is about > 
order of magnitude. The evolving spike gives the largest 
enhancement to the overall flux at the lower energy re- 
gion, while there is little enhancement at high energies 
close to m x . This is expected, as the spikes are most 
prominent just after formation at zbh , and gamma-rays 
emitted then have been redshifted to lower energies. Ex- 
cept for massive halos, most spikes today have died away 
and contribute very little to the gamma-ray flux. This 
also implies that the annihilation signal from the GC is 
not expected to vary significantly from the case of profiles 
without spikes. 

We also consider a case where gamma-rays are emit- 
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E (GeV) 

FIG. 5: Gamma-ray flux due to annihilation of neutralinos 
into two photons. Halos have mass between 10 5 - 10 14 M© 
and are distributed according to the Press-Schechter formal- 
ism. Flux from halo only are shown in black dotted line, and 
flux from halo and evolving spike are shown in red solid line. 
Spikes give the greatest contribution at low energies, i.e. at 
high redshifts. 



ted from annihilation of neutralinos into two photons. 
The photon flux is described as a delta function where 
dN(E)/dE = b J7 S(E - m x ), with b- n = 0.003 fl. Con- 
sidering only a delta function as the flux source is helpful 
to understand the enhancement due to the presence of 
spikes. The spike is expected to give the largest con- 
tribution around zbh, which today will be observed as 
gamma-rays with energy lower than m x . Figure [5] shows 
the flux from annihilation into two photons for halo only 
and halo and spike contributions, and indeed, the largest 
enhancement comes at low energies. We have again as- 
sumed that halos form at z = 18 and spikes form at 
zbh = 2 and used the M — Msmbh relation Eq. (JTJ)-(a). 
The steep enhancement for the spike's flux at E w 30 
GeV is due to our assumption of having a fixed SMBH 
formation epoch (zbh) an d only using the delta function 
for the gamma ray flux. 

To compare with existing literature, we introduce here 
a "boost factor" b max defined as in Ref. i.e. & max = 
min^EGRET^)/*^) Where $egret is the EGRET 
flux measurement (Ref. [28| for the diffuse background, 
and Refs. 29, 30] for the GC) relative to the energy bin 
Ei . We show in Fig. [B] the required boost factor for 3 
different cases: halo only (top solid line), and halo+spikc 
with 2 different assumptions for the Msmbh-M rela- 
tion (lower solid lines). The left panel shows a constant 
(crv) = 10~ 26 cm 3 s _1 with varying m x and the right panel 



shows a constant m x with varying (crv). The boost fac- 
tor for the GC is also shown as a comparision in dashed 
lines, where the HESS GC observation 3l| has also been 
considered. As one can see, for the Msmbh-M relation 
of Eq. lU)-(c), the required boost factor for the gamma- 
ray background is smaller than for the GC for most cases. 
We recall here that the spike contribution scales differ- 
ently with the particle physics parameters m x and (av) 
with respect to the halo only case, due to the satura- 
tion effects produced by annihilation itself. In order for 
annihilations to contribute significantly to the observed 
gamma-ray background, a boost factor of at least 2 orders 
of magnitude is thus required. This could in principle be 
achieved by steepening the halo slope in the innermost 
regions, for e. g. d ue to adiabatic compression of baryons 
(see e.g. Ref. [17| and references therein), or to the pres- 
ence of mini-spikes around intermediate mass black holes 
[l8l |32] ] . One should however bear in mind that astro- 
physical sources are expected to provide a significant, 
possibly dominant, contribution to the background. Fur- 
thermore, the estimate of the background measured by 
EGRET has actually been recently questioned by several 
authors. We discuss in the next section the uncertain- 
ties on the EGRET measurements and on the possible 
astrophysical intepretation, and in light of these uncer- 
tainties, we do not attempt to fit the background with 
a combination of particle physics and halo models, and 
limit ourselves to point out the importance of the role 
played by spikes in the estimates of the DM annihilation 
contribution to the extra-galactic flux. 

In Figs. HUH we have assumed a common redshift of for- 
mation for all SMBHs. We show in Fig. the dependence 
of the gamma-ray background on the parameter zb ti- 
the left panel shows the evolution of ((z) (Eqn. (|2"P|) ) for 
different values of zbh, and the right panel shows the 
gamma-ray background. Younger spikes give a greater 
contribution to the gamma-ray background because of a 
larger £(z). The normalization of the annihilation signal 
has a slight dependence on zbh, where small zbh values 
give larger contribution to the gamma-ray flux. This is 
expected as spikes that formed in earlier epochs evolve 
away with time. 

On the other hand, dependence on halo formation red- 
shift is negligible; changing the upper limit on redshift fro 
Eq. (fTT)) from 18 to 12 brings negligible change for both 
halo and halo+spike gamma-ray flux. The flux has a 
small dependence on the maximum halo mass: a 10 13 M 
limit lower the flux by < 20 %. Varying the minimum 
halo mass brings negligible change. 

All calculations so far assumed that spikes never expe- 
rienced a major merger, which could in principle signifi- 
cantly lower the DM density due to the souring effect of 
binary BHs (33[. To model the effect of galaxy mergers 
on the annihilation signal, we used the merger tree model 
of Ref. [2|| (Eq. (|A18|) ) , and assume that a galaxy merger 
occured and its spike destroyed at z m when its halo mass 
at zbh doubles, i.e., M(z m ) = 2M(zbh)- Figure [5] shows 
the effect merger has on the gamma-ray signals produced 
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(GeV) (<rv) (cm 3 s" 1 ) 

FIG. 6: Required boost factor to match the EGRET background measurements with DM annihilations, for the halo only case 
(top blue solid line); halo+spike with Msmbh — M relation (a) (middle red solid line); and halo+spike with Msmbh ~ M 
relation (c) (lower red solid line). Halos have mass between 10 5 - 10 14 M and are distributed according to the Press-Schechter 
formalism. Left panel: m x varies while the cross section is fixed at (av) = 10 -26 cm 3 s -1 . Right panel: (av) varies with a 
constant m x — 100 GeV. For comparison, we show the boost factor relative to the gamma-ray source at the GC (dashed line), 
using a NFW profile. Both EGRET and HESS observations are used. 
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FIG. 7: Left panel: £(z) for spikes (solid lines), assuming from bottom to top, zbh = 6 (black), 4 (green), 2 (red). The halo 
( z corresponding to the halo-only contribution is shown for comparion (dotted line). Right panel: Sensitivity of the predicted 
gamma-ray background to the formation redshift of SMBHs, zbh, compared with zbh = 2 case used in previous figures. 
The lower dotted line represents the contribution of the halo only, the four upper curves dashed and solid lines indicate the 
halo+spike contribution for the Msmbh-M relation (a) with, zbh = 6 (black), 4 (green), 2 (red solid), 1 (blue) from bottom 
to top. Halos have mass between 10 5 - 10 14 M Q and are distributed according to the Press-Schechter formalism. 
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FIG. 8: Gamma-ray background produced by DM annihila- 
tions in spikes including (dashed lines) and neglecting (solid 
lines) the effect of mergers. Three redshifts of formation have 
been considered for SMBHs; zbh = 2 (red, top lines), 4 
(green, middle lines), 6 (black, lower lines). The dotted line 
represents the halo contribution. Halos have mass between 
10 5 - 10 14 Mq and are distributed according to the Press- 
Schechter formalism. 



by spikes. The solid lines are the spike+halo contribution 
without mergers, and the dashed lines are the spike+halo 
contribution with merger taken into account. The dotted 
line is the halo contribution only, shown for comparison. 
Three redshifts are considered, zbh = 2, 4, 6, which gives 
a z m — 0.88, 2.12, 3.37, respectively. The reason for such 
small effect of mergers can be seen from the left panel 
of Fig \7\ most of the contribution from the spikes come 
right after its formation. By the time of the merger, 
most spikes would be already quite small and not give 
significant contribution to the gamma-ray background. 



IV. DISCUSSION AND CONCLUSIONS 

Different strategies have been proposed in the liter- 
ature to search for DM annihilation radiation. One 
of the most popular targets of indirect DM searches is 
the GC. The prospects for detecting gamma-rays from 
DM annihilations at the GC have been discussed ex- 
tensively in Refs. [H \M, H US [H, [M IH for DM 
cusps in the framework of different DM candidates, and 
in Refs. [H, [U [4l| for the case of a DM spike at the 
GC. An updated discussion in light of the recent discov- 
ery of a point source coincident with the GC, extending to 
very high energies can be found in Refs. [42|, 0, HH |45| . 



Current data do not allow a convincing interpretation 
of the gamma-ray emission as due to DM annihilation, 
while the properties of the gamma-ray emission appear 
consistent with those expected for an ordinary astrophys- 
ical source. However, although the DM intepretation of 
the gamma-ray source at the GC appears problematic, 
it can certainly be used in a conservative way to set up- 
per limits on the annihilation signal. Alternatively, one 
could search for the contribution of DM annihilations to 
the cosmological gamma-ray background, as discussed in 
Refs. 0, IS HIM El El. The "smoking-gun" in this 
case may come from the peculiar angular power spectrum 
predicted for this signal [-491 ]. 

The predicted gamma-ray background is usually 
compared with the extra-galactic emission measured 
by EGRET [5(|. The most convincing interpreta- 
tion in terms of conventional astrophysical sources in- 
vokes a large contribution from unresolved blazars (e.g. 
Ref. [5l|]), although this conclusion has been challenged 
by other authors (e.g. Refs. [H, HH). An additional 
contribution may arise from Inverse Compton scattering 
of electrons accelerated at shocks during structure for- 
mation [54], [H, , but this process can hardly account 
for the bulk of the background The EGRET extra- 
galactic background should however be treated with cau- 
tion, since it has been inferred (not measured) by sub- 
stracting the estimated Galactic foreground from the 
high latitude EGRET measurements. In a recent re- 
analysis of the EGRET data, Kehet et al. [HI], noticed 
that the high latitude profile of the gamma-ray data ex- 
hibits strong Galactic features and claimed that it is well 
fit by a simple Galactic model, obtaining an upper limit 
on the extra-galactic background 3 times stronger than 
previously assumed, and evidence for a much lower flux. 
In light of the large uncertainties associated with the data 
and with the contribution of conventional astrophysical 
sources, we conservatively consider the EGRET estimate 
as an upper limit to the actual gamma-ray background, 
and do not attempt to fit the data with an ad hoc com- 
bination of particle physics and halo properties. 

A comparison of the two strategies (GC vs. extra- 
galactic background) has been performed in Ref. 0], 
where it was shown that for ordinary cusps, and in par- 
ticular for an NFW profile, the prospects for detect- 
ing gamma-rays from the GC are always more promis- 
ing than for the gamma-ray background. Here we have 
shown that the situation changes, if we take into account 
the formation and evolution of DM spikes, which form 
due to adiabatic growth of SMBHs at the centers of DM 
halos. In fact, in this picture a spike inevitably devel- 
ops also at the center of the Galaxy, but it is rapidly de- 
stroyed by the combined effect of DM scattering off stars, 
and DM annihilations themselves. The enhancement of 
the annihilation signal is thus negligible [H, [TtJ ■ 

We have shown here that the opposite is true for the 
gamma-ray background. In fact, although all spikes are 
affected by the very same processes, the signal in this 
case receives contributions also from halos at high red- 
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shift, at a time when astrophysical and particle physics 
processes did not yet have the time to affect the DM 
density. As a consequence, the gamma-ray background 
from DM annihilations receives a substantial boost, so 
that its detectability is in some scenarios more promising 
than the case of a gamma-ray source at the GC. An addi- 
tional reason to consider the gamma-ray background as a 
valid alternative to GC searches, is that it is sensitive to 
the average properties of halos, whereas in the GC case 
one has to deal with a single realization that, as far as 
we know, may differ significantly from the average, given 
the significant scatter in the properties of halos observed 
in numerical simulations, and given the unknown history 
of the baryons. 

Several effects could further boost the annihilation 
background. One possibility is that DM halos undergo 
"adiabatic contraction" under the influence of baryons, 
thus steepening the original DM profile (see e.g. 17| and 
references therein), a circumstance that would lead to 
a similar boost of both the background and GC fluxes. 
Conversely, if mini-spikes of DM around intermediate 
mass black holes exist [HI, [32[ , this would have dramatic 
implications for the predicted annihilation background, 
while leaving practically unchanged the predictions for 



the GC M 
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APPENDIX A: PROPERTIES OF DM HALOS 

We present here our prescription to calculate DM den- 
sity profiles Ph(r), for halos of mass M and redshift z, 
and the dimensionless flux multiplier C,(z), that we used 
for the calculation of the gamma-ray background. For 
the sake of completeness, we explicitly write all the in- 
gredients of the calculation, with relevant references, also 
to allow the comparison with exisiting literature. 



1. Halo density profile 

The virial mass M of a DM halo can be expressed in 
terms of the virial radius 



M 



4tt 



■ Ac p c , 



(Al) 



where p c is the critical density and A c is the virial over- 
density, that for a flat ACDM cosmology can be approx- 
imated by [60j 

A c = 18tt 2 + 82{Q m (z) - 1) - 39(O m (z) - l) 2 . (A2) 

Here, £l m (z) is the matter density in units of p c . 

N-body simulations suggest that the density of DM 
follows a universal profile, usually parametrized as [l|, 



Ph(r) = po 



a + 0j 
'0 



(A3) 



where po and ro set the normalization of the density 
and radius, respectively. We adopt here the so-called 
Navarro, Frenk and White (NFW) profile [6l|, which can 
be obtained from the above parametrization with the fol- 
lowing choice of parameters a = /3 = 1, 7 = 2. In this 
case, the profile reduces to 



Pnfw{x) 



Po 



x (1 + x) 2 



(A4) 



where x = r/r$. The reference density po can be ex- 
pressed in this case as 



1 M 
Pa = TZ 3" 



1 



4?r r\ ln(l + c) - c/(l + c) 
where c is an adimensional concentration parameter 

r v 

c — — 



(A5) 



(A6) 



Although subsequent studies have refined the description 
of DM halos with respect to the NFW profile, and more 
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recent parametrizations of the innermost regions of galac- 
tic halos appear better suited to ca ptu re the behavior in 
the small-radii limit (see e.g. [62l . l63l|). we present our 
results for the NFW profile, which has emerged over the 
years as a benchmark model for DM annihilation studies. 
The scenario and prescriptions described here, however, 
can easily be extended to any DM profile. 

A convenient expression for c at redshift z for halo mass 
M{z) was derived in Ref. [64| from a statistical sample 
of high-resolution N-body simulations, containing ss 5000 
halos in the range 10 11 - 10 14 M Q : 



M{z) 



1 + z \M+(z = 0) 



-0.13 



(A7) 



The collapse mass M*(z) defines the mass that collapses 
to form at halo at redshift z. For ACDM cosmology, 
M*(z = 0) ~ 1.5 x lO^h' 1 M Q . A detailed way of 
obtaining M* is described in Sec. IA 31 Concentration pa- 
rameters of halos less massive than ~ 10 10 M sun have not 
been robustly measured in N-body simulations, due to 
the required mass and force resolution. For concreteness, 
we assume Eq. (|A7[) to apply for any halos. The validity 
of extrapolating Eq. (|A7[) to small halos should be tested 
against future simulations. 



2. Halo mass distribution 

In this subsection, we discuss another key ingredient 
for the calculation of the cosmological annihilation flux. 
We need now to specify the mass function of DM halos, 
i.e. the number of objects of given mass M. The mass 
function of halos are expressed in a universal form [65[ 



dn _ pm_ dv_ . 
dM ~ M dM J ^ V ' 



(A8) 



where pM = Pc^m is the comoving matter (background) 
density and a new variable v is defined 



1.686 



D(z) a(M) 



(A9) 



The linear equivalent of the over density at collapse for 
spherical collapse model is 1.686, a(M) is the rms density 
fluctuation in a sphere with mass M, and D(z) is the 
linear density growth rate. 

In this way, f{v) for spherical collapse is expressed as 



(AlO) 



A better fit to number density of halos can be obtained 
with the ellipsoidal collapse model derived by Sheth and 
Tormen [66] ]; 

/m = (i + ^) {~ e ~ v,2/2 > ( An ) 



where v' = yfav. Numerical fits to simulations give a = 
0.707, p — 0.3 [67], and A is obtained by normalising 

/(")• 

a(M) is related to the power spectrum P(k) by 



a{Mf cx / d 3 kW(kR) 2 P(k). 



(A12) 



where k is the wave number and W(kR) is the filter func- 
tion with R being the radius enclosing mass M. We 
choose a top-hat function for W(kR); 



W(y) = o (sin y - y cos y) 

We adopt here a power spectrum 

P{k) cx k n T(k) 2 , 



(A13) 



(A14) 



where the Bardccn-Bond-Kaiser-Szalay (BBKS) transfer 
function is used for T(fc): 



ln(l + 2.34o) , 
T{k) = ——— [1 + 3.89q- 



2.34q 
(m.lq) 2 + (5.46g) 3 



(6.71 g ) 



41-1/4 



(A15) 



with q = k/(Q, m h 2 ) (Mpc). P(k) and a(M) are normal- 
ized by simulating a(M) in spheres of R = 8/h Mpc, 
commonly known as erg. We use ag — 0.9 from the con- 
cordance model. 

Finally, we provide an expression for the linear growth 
rate D(z) [69(, which can be expressed as 



D(z) 



i g(z) 



1 + z g(0) 



(A16) 



where 



g(z) = - tt m (z) x fi m (z) 4/7 - 



70 



(A17) 



3. Evolution of halos 

In order to assign the appropriate BH mass to a host 
halo, we need to evolve back in time the halo mass, and 
calculate its mass at the redshift of formation of the 
SMBH. We follow the semi-analytic study of halo evo- 
lution with merger trees, carried out in Ref. [26j , and 
express the halo mass at zq as a function of an earlier 
redshift Z\ 



M(zq) = M(zi) exp 



S 



1 



1 + z c \ 1 + z 



- 1 



(A18) 



The free paramter S is proportional to the logarithmic 
slope of accretion rate, dlogM/da where a = 1/(1 + z), 
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and is usually set to 2. The collapse redshift z c is defined We solve the above equation to find z c for each M. 
according to Ref. where the collapse mass is a By setting v = 1, we have 
fixed fraction of the halo mass M 

1.686 _ 1.686 

M*(z c ) = fM. (A19) ° {Zc) ~ a(M*(z c )) ~ tr(O.OlM) ' ( } 



In ACDM cosmology, / is typically 0.01. M+(z) is defined Hence z c is obtained by interpolating the expression 
such that a(M ir (z), z) = 1.686, or equivalently, v = 1. -D(z c ) of Eq. (|A16jl . 



